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1. Introduction 

The Weierstrass criterion is one of the most elegant and powerful tools of classical 
mechanics, in that it allows qualitative analysis of one-dimensional conservative 
motions, once the zeros of a suitable function, the Weierstrass function, are de- 
termined. The usefulness of this criterion also lies in the fact that it can also be 
successfully generalized to many problems outside of classical mechanics, see e.g. 
the recent papers 1 112131 . 

In particular, the focus of this paper is to study the qualitative behavior of 
Lemaitre- Tolman-Bondi (LTB) models, endowed with a non-null cosmological con- 
stant A, by assuming a Weierstrass approach. @ 

a A similar approach has been developed in II 121 . but in that case only LTB models with null A 
were considered. 



2 /. Bochicchio, S. Capozziello and E. Laserra 



The importance of such Universes in relativistic cosmology is mainly based on the 
following observation: even if the Friedmann-Robertson- Walker (FRW) modelsE) 
are widely accepted to describe the cosmic evolution, there are several objections to 
such models which include singularities, horizons, observed inhomogeneitics start- 
ing from galactic scales up to galaxy supcrclustcrs (e.g. Virgo superclustcr) [4]. 
In order to overcome some of the difficulties faced by the FRW models, Lcmaitre 
[5], Tolman [B], Bondi [7] and others have considered inhomogeneous spherically 
symmetric dust models. Hence, the main aim of the LTB approach is to encom- 
pass cosmic inhomogeneities, at both large and small scales, with the overall cosmic 
dynamics JHJ. Besides, being the simplest inhomogeneous solutions of the Einstein 
equations, it is relatively easy to work with them. Furthermore, the interest in such 
models is recently increased due to the fact that some of them can be designed to 
satisfy several observational requirements |9|10|llj . 

As such, through the Weierstrass method, we are able to qualitatively describe 
and classify the possible kinds of evolution of the LTB-models. In particular, in the 
peculiar case of the FRW-models, Fricdmann discussion [12|13] can be reduced to 
a straightforward application of the Weierstrass method. 

The paper is organized as follows: in Sec. 2, the main features of LTB Universes 
are reviewed. A discussion on the qualitative study of the evolution of the r-shells 
in General Relativity through the Weierstrass method is given in Sec. 3. Section 4 
is devoted to concluding remarks. 

2. The Lemaitre Tolmann Bondi Models 

Let us give now a brief summary of the main features of LTB-models according to 
[1417115116] . We will consider a dust system C which, during its evolution, generates 
a Ricmannian manifold, with locally spatial spherical symmetry around a physical 
point O;0 the metric can then be given the Levi Civita's form |17j : 

ds 2 = A 2 {t,r)dr 2 + R 2 (t,r) (d 9 2 + sin 2 ddip 2 ) - c 2 dt 2 , (1) 

where t is the proper time of each particle and r, (9, ip are co-moving Levi-Civita's 
curvature spherical coordinates: we can interpret R(t, r) as the intrinsic radius of 

the O-sphere S{r) at time t so that — ^ represents, at any point, the Gaussian 
curvature of the geodesic sphere with its centre at O. 

We consider now the initial space-like hypersurface V3 (with equation t = 0) 
and call r-shells the set of particles with co-moving radius r (i.e. the dust initially 
distributed on the surface of the geodesic sphere with center at O and radius r 
(O-sphere) S(r)). According to [14115116] . we assign, at each particle of a r-shell, 
the initial intrinsic radius r = R(r, 0) as the radial co-moving coordinate. Hence 

b That is the homogeneous spherically symmetric dust models. 

c See [17] for a precise definition of locally spatial spherical symmetry around a physical point O. 
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the metric of the initial spatial manifold V3 takes the form 

dcr 2 = a 2 (r)dr 2 +r 2 {dd 2 +sin 2 6»d^ 2 ) , 



(2) 



where a(r) = A(r, 0). 

Now let us take into account the gravitational field equations with the cosmo- 
logical constant A and the conservation equations: 

8ttGat 

1 1 a0 1 

(3) 

V a T| = , 



Gap + A g a £ 



where Gn is the Newton gravitational constant. 

It has been shown [18] that it is possible to break the corresponding Cauchy problem 
into two separate intrinsically formulated invariant problems: the problem of initial 
conditions and the restricted problem of evolution. By taking into account [14171 
I18I15I16] . we can immediately translate the restricted evolution problem into the 
following equations^ 



' A{t,r) = a(r)R'(t,r) 
R 2 

fj,(t,r) 



a 2 (r) 



2G N M(r) A 2 
R 3 



(4) 



R'(t,r)R 2 (t,r) 

where a dot denotes differentiation with respect to t and a prime differentiation with 
respect^ to r, n(t,r) is the mass density (fio — /i(r, 0)) and M(r) is the "Euclidean 
mass r 



'0 



M(r) = 4?r / fi (s)s z ds . 
Jo 



(5) 



If the initial mass density is constant, equation ((4]) 2 reduces to the Friedman 
equation and M (r) is constant. 

If, analogously to [T^] , we introduce the function 



e(r) 



a 2 (r) 



(6) 



which represents the percentage deviation of a 2 (r) from the Euclidean value a 2 = 1, 
Eq. (QI2 becomes 



i? 2 



e{r) c 



2 Gat M(r) 



A 



R 



-R 2 



(7) 



d These equations differ from the previous 1 151 161 only for the additional term ^ B? . 

6 In that M(r) would be the mass of dust contained within S(r), if the initial O-sphere V3 was 

Euclidean. 
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Remark 1. The function e(r) and the initial metric coefficient a 2 (r) are connected 
to the initial data by the equation 



which gives an initial constraint. 

3. Generalization of the Weierstrass Criterion to the LTB Models 
3.1. The Weierstrass equation 

We will use a generalized form of the Weierstrass criterion,0 which can be applied 
to many problems outside of classical mechanics too. Let's consider the first order 
differential equation 

x 2 = $(» , (9) 

that we may call Weierstrass equation with Weierstrass function <fr(x) (see [20] ) . 
Equations of the form ^ arc frequently encountered in classical mechanics. For 
example, the natural motions of a material point of mass m, subjected to a conser- 
vative force deriving from the potential energy V(x) and having mechanical energy 
E, are described by the Weierstrass equation 

x 2 = -[E-V(x)] , (10) 
m 

whith the Weierstrass function 

®(x) = — \E-V(x)] . 
m 

In addition the Weierstrass Eq. ([9]) translates into the double equation 

^=±V¥£), (11) 
which can be integrated by separating the variables 

t(x)=± f -^=+*o (12) 

where wc choose the sign ± in agreement with the sign of the initial rate Xq, 

±1 = <f>(x ) . 

The importance of the Weierstrass approach is mainly based on the fact that 
it is possible to obtain the qualitative behavior of the solutions of a Weierstrass 
equation, without integrating it. Precisely, the zeros of the Weierstrass function 
have a leading role, in fact the solutions of the Weierstrass equation are confined 
in those regions of the x-axis where x 2 > 0, hence the Weierstrass condition 



> (13) 



f For a detailed exposition of the Weierstrass criterion see e.g. [201 § 1.3 p. 23]. 
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must be fulfilled: the solutions of the Weierstrass equation are confined in those 
regions where the Weierstrass function is non negative. 

These regions are unlimited or are limited by the extrema of the definition interval 
(eventually +00 or —00) and by the eventual zeros of the Weierstrass function, 
which are called barriers, because they cannot be crossed by the solution x(t), so 
they split the range of possible values for x into allowed and prohibited intervals. 
Now let's consider the zeros of $. 

If a barrier aig is a simple zero of that is such that 

then it is called an inversion point Xi, because the motion reverses its course after 
reaching it (see e.g. |2QI3j ). 

If a barrier is a multiple zero, that is such that 

$(xb)=0, &(x b )=0, 

then it separates two allowed intervals and it is called a soft barrier xs [3]. 

A soft barrier is also called an asymptotic point, because it takes an infinite time 

to reach it. In fact at a soft barrier the integral (TT2")) diverges (see e.g. |20I3| ). 

Finally we recall that an asymptotic point is also an equilibrium point (see e.g. 

[20]). 

So, once these zeros are found, the qualitative behavior of the solutions of the 
Weierstrass equation ((9J is completely determined. 



3.2. The Weierstrass criterion for the evolving r-shells 

In order to qualitatively study the behavior of r-shells in the case of non-null 
cosmological constant, let us analyze Eq. (j4|)2 which determines the evolution of 
the material continuum through the Weierstrass method. 

Now we will focus our attention on a given single r-shell (that is we will consider 
r as a fixed parameter), so we can regard the intrinsic radius R as a function of 
the time only, x(t) = R{t;r), and e(r) ,M(r) as constant; then Eq. ([7]) becomes a 
quadratic differential equation 

x 2 =ec 2 + + -x , (14) 

x 3 

that is a Weierstrass equation with Weierstrass function 

, . 9 2 G N M(r) A 9 , . 

$(x;r =erc 2 + ^ + - x 2 , (15 

x 3 

depending on the parameter r. Eq. (|14[) translates into two equations, depending 
on the parameter r, 

f t = ±VH^) - ± ] /e(r) * + , (16) 

g Once solved Eq. $4^2, Eqs. 1 5 (HJ3 can be immediately solved. 
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where we have to choose the sign ± in agreement with the sign of the initial rate 

•2 ^, x / s. 2 2 Gn M(r) A 2 
x 2 = $(a; ; r) = e(r) c 2 + — ^ + -r 2 , 



Remark 2. In consequence of Poincare's theorem on the analytic dependence of a 
solution on a parameter, the solutions of each of the Eqs. (|T51) depend analytically 
on the parameter r in each point where depends analytically on x and r (see 
e.g. [U). 



The zeros of the Weierstrass function (|T5j) depend not only on the sign of A, but 
obviously on the sign of e too, so the evolution of the r-shells will depend on both 
signs, the sign of e and the sign of A. These different situations will be analyzed 
in details in the following sections, where the zeros of the Weierstrass function are 
obtained by finding the positive real roots of the third degree equation 

Aa; 3 + 3ec 2 a; + 6GjvA/ = . (17) 

To highlight the role of the cosmological constant, in the following we will consider 
Eq. (|T3|) in the form: 

x 2 = l -x 2 \k + W{x)] , (18) 



. , 3ec 2 6G N M 
W(x) = + . (19) 



where 



Then the Weierstrass condition (|13|) holds when 

A > -W(x). (20) 



Remark 3. Once W(x) is introduced, a barrier xb 0) is a soft barrier xs iff 
W'(xb) = 0. 

In fact from Eqs.© and (fT8]) it follows that 

= \x [A + W(x)] + \x 2 W'{x). 

Supposing $(.t b ) = then A + W{x B ) = 0; so, since $'(cc B ) = \x 2 W'(xb), 
&'(xb) = => W'(xb) = and viceversa. 

We remark that, by considering the Weierstrass function <J>(x) in the peculiar 
case of a FRW-Universe, it is easy to re-obtain the Friedmann discussion [12] . 
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3.2.1. Null cosmological constant A 

For sake of completeness, in this section we briefly recall the qualitative behavior 
of LTB-models with A = 0, since they have already been studied in |1I2) . 
For negative values of e, there is only one barrier 

2GM 

XI = J" ' 

which is a simple zero of the Weierstrass function, that is an inversion point. So, 
if a r— shell is initially expanding, it will go on expanding until the intrinsic radius 
reaches the value x/; then it will contract back from xi towards the center of 
symmetry O until it collapses in a finite time (see [1]). 

When e is null (Euclidean case) the inversion point goes to infinity, so if a i — shell 
is initially expanding, it will go on expanding without limit, approaching the null 
expansion rate (see [1]). 

Finally, a similar result can be achieved for positive e: there are no barriers, 
so if the r-shell is initially expanding it will go on expanding with decreasing rate 
approaching the limit value x\ = V ec 2 (see [T]). 




Fig. 1. The Weierstrass function <J>(x) for A = in three different cases. More precisely 

i) when e < 0, the function has a simple zero at Xj (i.e. there is one inversion point Xj) and 
the motion V < X < x r is possible; 

ii) when e = 0, the inversion point goes to infinity and the motion is possible V x > 0; 

iii) when e > 0, there are no barriers: the motion \/ x > is possible. In this case the Weierstrass 
function admits the red line &(x) = ec 2 as asymptote. 
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3.2.2. Positive cosmological constant A 

The positivity of the cosmological constant implies an open model for null and 
positive values s(r), while a different and more complex situation is obtained for 
negative values of e(r). 



A is positive and e is positive or null 

The evolution of LTB-models in these two cases is very similar. 
In fact when e = (Euclidean case), the function (fT9)) becomes 

w , , 6G N M 

so, from (HOD, if A > then $(x) > Vz > 0. 

Instead when e > 0, the function (fTi?]) becomes 

„,. , 3ec 2 6G N M n 

w ( x ) = — + —^r- > 

and, from ((20]), if A > then $(ar) > Vx > 0. 

Hence, in both cases, there are no barriers and we have a monotonic expansion: if 
the r -shell is initially expanding, it will go on expanding without limit, otherwise it 
will collapse. 




Fig. 2. The Weierstrass function <E>(a;) for positive A. Cases e = and e > are analyzed. Here 
there are not inversion points: the motion Vi > is possible. 
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Positive A and negative e 

When e is negative, the function (fT9]) becomes 



+ (21) 



and it has a minimum at0 



3G N M 

xm = i i n ■ (22) 



Moreover, from Remark [31 it is clear that when A assumes a suitable critical value 
A c defined as 

A c = -W( XM ) = Jfj^L , (23) 



the Weierstrass function (|15[) has a soft barrier in the point Xm, that we can also 
write as 



X M 



3G N M 

~^7~ 



Hence when A assumes its critical value A c , the evolution of r-shells in LTB-models 
is static and stable. 

Precisely the behavior of each r-shell depends by the initial conditions: their evo- 
lution is really static if the initial intrinsic radius equals xm', on the other hand, 
the static situation is a limiting situation: if xo ^ Xm, the r-shcll go on expanding 
asymptotically approaching to the static model with xm as intrinsic radius. 

If A > A c , there are no barriers, hence if the r-shell is initially expanding 
(xo > 0), it will go on expanding without limit. 

Finally, when < A < A c , the Weierstrass function (|T5j) admits two simple 
zeros 



,Jc 2 a + 47T /He 2 a 

x\ = 2\ — ; — cos and x-i = 2\ — - — cos — , 

V A 3 V A 3 

where a is defined by 

tana = , 



-Q 

W i th q = and A = 9G* N M>A-\ e f c ° J 

So when the initial intrinsic radius is less that x\ and the r-shell is initially expand- 
ing, it will go on expanding until the intrinsic radius reaches the maximal expansion 

11 In fact its first a: -derivative = ~^~jr 1sg n m j g nu [[ wnen x = xm- In addition, 

2j £ |5 c 10 

evaluating the second derivative of W(x) in m, we obtain the positive value 4 4 . Hence 
xm is a minimum point for the function W(x). 

'Note that A is negative since < A < A c . Moreover, since the tangent of a is negative, ct £ ( ^ , 7r) , 
hence xi < X2 ■ 
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point xi ; then it will contract back from x\ towards the center of symmetry O until 
it collapses in a finite time; on the other hand, when the initial intrinsic radius is 
grater than x~2 and if the r-shell is initially expanding, it will go on expanding. 
Note that the open or closed evolution of the r-shell strongly depends from the 
initial conditions. 

Moreover, a particular situation is observed when the initial intrinsic radius xo is 
such that X\ < xq < x~i , since in this case the Weierstrass function is negative: the 
initial intrinsic radius can't belong to the interval (xi,X2). 




Fig. 3. The Weierstrass function &(x) for positive A and negative e. In this case the behavior of 
the function depends on A. More precisely 

i) when A = A c , the function <E>(a;) has a multiple zero at x M . An equilibrium position or a limiting 
position corresponds to the soft barrier. In the first CelSGj X — X *r IS the unique possible position. 
In the second case, the motion V < x < x M is possible; 

ii) when < A < A c , the function <&(x) has two simple zeros at x\ and X2 (i.e. there are two 
inversion points x\ and X2). The motion is possible V < x < x± and Vx > X2\ 

iii) when A > A c the Weierstrass function has no zeros. The motion Vi > is possible. 



3.2.3. Negative cosmological constant 

Let us focus now on the qualitative behavior of r-shells for negative value of A. 
In this case, the zeros of the Weierstrass function may only be simple zeros, which 
correspond to inversion points. Even if the expression of the maximum value reached 
by the intrinsic radius depends by the sign of e, for each analyzed situation (e = 0, 
e < and e > 0) the model is closed. 
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Table 1. Evolution of of LTB r-shells for positive cosmological constant and different 
signs of e(r). 



Sign of e(r) 



Evolution of r— shells 



A > 



€ = 



e < A 



e < 0, A > A c 



e<0,0<A<A c 



e > 



Open 



Static solution as effective or limiting solution 
In the second case the evolution is open 



Open 



Closed if xq < x\ — 2\l 
Open if xq > x~2 — 2 



Open 



Tabic 2. Evolution of LTB r-shells for negative cosmological constant and different signs of 
e(r). 



Sign of e(r) 



Evolution of r-shclls 



A < 



e = 



e < 



Closed 



Closed 



6 > Closed 

(the maximum value of the intrinsic radius depends also by A) 



Negative A and null e 

When e = (Euclidean case), we find the barrier 

_ ( 6G N M \i 
Xs ~ { -A ) ' 

which is a simple zero of the Weierstrass function (|15|) and corresponds to an 
inversion point. We have $(x) < for x > x si hence the evolution is possible 
for all values < x < x s . 

When e = and A < 0, if the r-shell is initially expanding, it will go on expanding 
until the intrinsic radius reaches the invrsion point x s ; then it will contract back 
from x s toward the centre of symmetry O until it collapses in a finite time. 
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XLimit Xi X 3 H I(l,2) 

Fig. 4. The Weierstrass function &(x) for negative A in three different cases. More precisely 

i) when e = 0, the function &(x) has a simple zero at x s (i.e. there is one inversion point x s ) and 
the motion V < x < x s is possible; 

ii) when e < 0, the function &(x) has a simple zero at x t (i.e. there is one inversion point x ; ) and 
the motion V < I < I, is possible; 

iii) when e > the Weierstrass function has a simple zero. The value of this zero depends on A 
as above specified. It can be individuated by xj l or by xj 2 (i.e. there is one inversion point xj 1 
or Xf 2 ) and the motion V < X < xj 1 (or clearly V < X < ^7 2 ) is possible; 

finally, when e < — 1 it is represented the limiting case e = — 1 when x^ imit is the lower bound 
for the set of inversion points x t . 



Negative A and negative e 

When A and e arc both negative, there is a barrier 



3G N M 




|e| 3 c 6 


|A| +' 


i iai 2 


|A|3 


3GivM 




|e| 3 c 6 


|A| 


V IAI 2 " 


|A|3 



+ 



(24) 



The barrier is a simple zero of the Weierstrass function, that is an inversion point , 
so the model is again closed: if the i — shell is initially expanding, it will continue to 
expand until the intrinsic radius reaches its maximum value x t ; then it will contract 
back from x I toward the center of symmetry O until it collapses in a finite time. 
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Negative A and positive e 

Let's put, for sake of commodity, 



9G 2 N M 2 



|A| 2 |A|* 

When 



|A|> 



£ 3 C 6 



9G N 2 M 2 ' 

we have A > and there is a barrier 

(3G H M i\* /3C„M ., 

*"-(n3r +A, J + \—- A 

which is a simple zero of the Weierstrass function, so corresponding to an inversion 
point. 

Instead, when 

£ 3 C 6 

|A|< 9G^' 
A is strictly negative; let's put, for sake of commodity, 

^A>r 



a = arctan 
then there is a barrier 



3G N M 



„ / 3ec 2 a 

^ =2 V^vT C0S 3 ' 

which is again a simple zero of the Weierstrass function, so corresponding to an 
inversion point. 

Hence, if the r -shell is initially expanding, it will continue to expand until the 
intrinsic radius reaches the value xi 1 or xi 2 ; then it will contract back from xi 1 or 
xi 2 toward the center of symmetry O until it collapses in a finite time. 



4. Concluding Remarks 

In this paper, we have analyzed a cosmological scenario based on the spherically 
symmetric dust solutions of the Einstein equations, that is LTB models with non- 
null cosmological constant, through an approach strictly related to the Weierstrass 
method of Classical Mechanics. In our case, it allows a systematic analysis of LTB 
models according to the different signs and values of e(r) and A. 

In particular, while we can easily classify the LTB evolution with null cosmolog- 
ical constant |l|2j . namely the Universe is open when e > and closed if e < 0, the 
situation in presence of cosmological constant is much more complex. First of all, 
one has to study the evolution for positive and negative values of A and, for each of 
these, relate the evolution to the sign of e. In particular, when A > 0, for e > 0, the 
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same global behavior of LTB models with null cosmological constant is achieved: 
the Universe is spatially open. In the other cases (e < 0) the Universe can be open 
or closed and this feature strongly depends on A and on the initial conditions. The 
various dynamical cases arc summarized in Table I. 

On the other hand, when A < 0, for any value of e, LTB Universes present the 
same evolution: they are always spatially closed. In particular, the maximum value 
reached by the intrinsic radius depends only on e when e < and also on A when 
e > 0. The dynamical behavior is more simple and summarized in Table II. 

In addition note that, in the cases corresponding to spatially closed Universes, 
we are not only able to select such a characterizing feature, but also to obtain the 
exact value of the maximum that the r-shell can reach. 

As concluding remark, it is worth noticing that the method outlined here can 
be, in principle, applied any time that cosmological dynamical system can be recast 
along the Weierstrass function and an effective cosmological constant can be defined. 
For example, it could be possible to apply this approach to cosmological models 
including scalar fields |22j . 
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